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Abstract 

Distributions  in  reactant  species  concentration  in  a  PEMFC  due  to  local  consumption  of  fuel  and  local  transport  of  water  through  the  membrane 
cause  distributions  in  current  density,  temperature,  and  water  concentration  in  three  dimensions  in  a  PEMFC.  These  distributions  can  lead  to 
flooding  or  drying  of  the  membrane  that  may  shorten  the  life  of  an  MEA.  Changing  the  cell’s  flow-field  pattern  to  distribute  the  gas  more  evenly 
is  one  method  of  minimizing  these  stresses.  This  paper  investigates  how  200  cm2  serpentine  flow-fields  with  different  number  of  gas  paths,  and 
thus  different  gas  path  lengths,  affect  performance  and  species  distribution.  The  results  show  how  the  local  temperature,  water  content,  and  current 
density  distributions  become  more  uniform  for  serpentine  flow-field  designs  with  shorter  path  lengths  or  larger  number  of  channels.  These  results 
may  be  used  to  develop  universal  heuristics  and  dimensionless  number  correlations  in  the  design  of  flow-fields  and  stacks. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

During  operation  of  a  proton  exchange  membrane  fuel  cell 
(PEMFC),  hydrogen  and  oxygen  are  consumed  as  the  fuel  and 
air  travel  from  inlet  to  outlet  leading  to  distributions  in  reaction 
rate.  This  distribution  of  the  reaction  causes  gradients  in  tem¬ 
perature  and  water  production  over  the  area  of  the  cell.  These 
secondary  effects  can  feedback  to  change  the  reaction  rate  and 
can  also  lead  to  ternary  effects  such  as  local  flooding  when  the 
local  partial  pressure  of  water  exceeds  the  saturation  pressure 
for  water  at  the  local  temperature.  Other  ternary  effects  include 
transport  distributions  of  water  across  the  membrane  between 
the  anode  and  cathode  due  to  electroosmosis  and  diffusion  as 
well  as  variation  of  membrane  conductivity  due  to  changes  in 
the  local  water  activity  on  the  anode. 

Many  researchers  have  studied  the  problem  of  water  manage¬ 
ment  inside  PEMFCs  for  both  steady  state  and  transient  oper¬ 
ation  [1].  However,  the  effect  of  flow-field  design  on  PEMFC 
performance  has  received  less  attention.  Fell  et  al.  [2]  used  a 
single-phase  isothermal  CFD  model  of  PEMFC  segments  to  test 
performance  of  experimental  flow-field  designs  under  differ¬ 
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ent  humidity  conditions.  Greenway  et  al.  [3]  used  a  two-phase, 
non-isothermal  model  to  study  10  cm2  serpentine  laboratory- 
scale  flow-field  patterns  with  different  flow  path  lengths  and 
inlet/outlet  configurations  to  study  cell  performance  and  distri¬ 
bution  uniformity.  They  also  reported  the  effect  of  gas  diffusion 
layer  (GDF)  properties  on  these  flow-fields.  Oosthuizen  et  al. 
[4]  investigated  mass  transport  under  the  ribs  in  a  PEMFC  using 
a  single  phase  non-isothermal  PEMFC  model  in  the  finite  ele¬ 
ment  code  FIDAP  to  study  channel  to  channel  transport  through 
the  GDF  for  different  GDF  permeability.  Experimentally,  Fi  et 
al.  [5]  compared  a  serpentine  pattern  to  a  matrix  of  individual 
graphite  squares  and  Yoon  et  al.  [6]  tested  different  rib  widths  on 
an  80  cm2  cell  to  see  the  effect  on  performance.  Commercial- size 
cells  (~ 200-800  cm2)  have  more  pronounced  effects  from  pres¬ 
sure  drop,  flooding  distribution,  and  non-uniformity  in  temper¬ 
ature  and  current  distributions  due  to  manifolding  constraints. 
These  effects  cause  stability,  durability,  and  performance  prob¬ 
lems  for  PEMFC  stacks.  Shimpalee  et  al.  [7]  explained  effects 
seen  in  a  patented  480  cm2  cell  under  different  cathode  humidity 
conditions,  but  did  not  investigate  different  flow-field  patterns. 

The  main  goals  of  flow-field  design  are  to  increase  uniformity 
of  the  current  and  temperature  distributions  at  the  operating  con¬ 
ditions  of  interest  while  maintaining  or  improving  polarization. 
Studies  of  different  flow  patterns  will  give  an  idea  of  how  to  opti¬ 
mize  the  flow-field  design  for  a  PEMFC  stack.  This  optimum 
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Fig.  1.  Schematic  of  the  typical  PEM  fuel  cell  used  for  testing  of  membrane  and  electrode  assembly  (MEA). 


depends  on  temperature  changes  inside  the  system  resulting 
from  the  electrochemical  reaction  and  water  phase  change  and 
it  has  implications  for  the  design  of  the  cooling/heating  system. 
Here  we  consider  the  flow-field  design,  temperature  change,  and 
effects  from  the  phase  change  of  water  on  the  performance  of  a 
commercial  size  (200  cm2  reaction  area)  PEMFC  performance 
at  steady  state  operation. 

Fig.  1  illustrates  a  schematic  of  the  typical  PEMFC.  The 
cell  consists  of  two  graphite  plates  with  machined  flow  chan¬ 
nels  separated  by  a  membrane/electrode  assembly  (MEA).  The 
MEA  is  a  sandwich  of  a  membrane  between  two  electrodes  with 
dispersed  platinum  catalysts.  In  this  paper,  we  will  focus  on 
the  five  different  200  cm2  flow-field  patterns  shown  in  Fig.  2. 
They  are  (a)  3 -channel  serpentine,  (b)  6-channel  serpentine,  (c) 
13-channel  serpentine,  (d)  26-channel  serpentine,  and  (e)  26- 


channel  complex  serpentine.  Geometries  (a)-(d)  have  similar 
patterns,  but  different  gas  path  lengths.  Geometry  (e)  has  been 
investigated  experimentally  by  Buchi  and  Neto  [8]  and  it  has  the 
same  path  length  as  the  26-channel  serpentine  flow-field.  This 
complex  geometry  has  been  included  into  this  study  to  exam¬ 
ine  the  effects  of  the  flow-field  pattern.  The  solutions  from  the 
model  will  be  shown  and  analyzed  for  these  different  flow-field 
configurations.  The  pattern  of  this  distributed  network  and  the 
effect  on  pressure  drop  will  be  studied  and  generic  design  prin¬ 
ciples  will  be  developed  that  we  hope  could  benefit  the  industry. 
The  goal  of  this  research  is  to  add  to  the  knowledge  base  to 
produce  generic  design  information  for  the  flow-fields  that  can 
be  applied  to  the  fuel  cell  stacks.  Heuristic  models  are  needed 
because  modern  day  computers  even  with  parallel  processing 
are  slow  to  compute  optional  stack  designs  without  some  intel- 


(C) 


Y 

L _ X 


(B) 


(D) 


(E) 


Fig.  2.  Flow  filed  patterns  of  anode  and  cathode  on  200  cm2  PEMFC:  (A)  3-channel  serpentine  flow-held;  (B)  6-channel  how-held;  (C)  13-channel  how  held;  (D) 
26-channel  how  held;  (E)  26-channel  complex  how-held. 


400 


S.  Shimpalee  et  al.  /  Journal  of  Power  Sources  160  (2006)  398-406 


ligence  in  the  set  of  flow-fields  considered.  This  study  is  a  part 
of  our  main  research  goal  of  developing  universal  heuristics  and 
dimensionless  number  correlations  that  can  be  used  to  design 
flow-fields  and  stacks.  It  is  assumed  that  the  development  of 
these  design  techniques  with  CFD  will  require  verification  (of 
a  subset  of  model  predictions)  but  that  proper  design  of  exper¬ 
iments  will  expedite  the  development  of  the  heuristics.  That  is, 
CFD  experiments  will  be  less  costly  than  build-and-test  exper¬ 
iments  and  they  can  be  correlated  to  yield  design  heuristics  for 
the  next  generation  of  PEMFCs. 

2.  Model  development 

The  equations  solved  in  this  work  included  the  conserva¬ 
tion  of  mass,  the  Navier-Stokes  equations,  the  species  transport 
equations,  the  energy  equation,  and  the  water  phase  change 
model.  Source  terms  for  these  equations  used  in  the  model  were 
taken  from  Ref.  [15].  In  the  water  phase  change  model,  when 
the  local  activity  of  water  exceeded  1 .0  water  vapor  was  con¬ 
densed  to  form  liquid  water  until  the  local  activity  equaled  1.0. 
Conversely,  if  liquid  water  was  present  and  the  local  activity  of 
water  dropped  below  1 .0  then  liquid  water  was  evaporated  until 
the  local  activity  equaled  1.0.  When  liquid  water  condensed  in 
the  region  adjacent  to  the  electrode  surface,  it  is  assumed  to  form 
a  liquid  film  on  the  electrode  surface.  In  the  regions  where  there 
was  a  liquid  film  on  the  electrode,  hydrogen  and  oxygen  were 
required  to  dissolve  in  the  liquid  film  and  diffuse  through  the 
film  to  the  electrode  surface  in  order  to  react.  Henry’s  law  was 
used  to  calculate  the  solubility  of  the  gases  in  the  liquid.  The 
thickness  of  this  liquid  film  depended  on  the  rate  of  condensa¬ 
tion/evaporation  and  the  production  of  water  by  electrochemical 
reaction  as  discussed  in  the  appendix  of  Lee  et  al.  [9]. 

A  control  volume  technique  based  on  a  commercial  flow 
solver,  STAR-CD  3.150A026,  was  used  to  solve  the  coupled 
governing  equations.  This  software  was  used  with  an  add-on 
tool  called  expert  system  for  PEMFCs  (ES-PEMFC)  that  pro¬ 
vided  the  source  terms  for  the  species  transport  equations,  the 
phase  change  equations  for  water,  and  the  heat  generation  equa¬ 
tions  [13].  Also,  ES-PEMFC  accounted  for  the  flux  of  protons 
and  water  across  the  membrane  [13]. 


The  PEMFCs  simulated  in  this  work  (shown  in  Fig.  2)  con¬ 
sisted  of  two  flow-field  patterns  (upper  is  anode  and  lower  is 
cathode)  separated  by  GDLs  and  an  MEA.  Table  1  summarizes 
the  geometry  details  for  each  of  the  models  and  properties  used 
in  the  simulation.  The  length  of  the  flow  path  for  each  geome¬ 
try  was  inversely  proportional  to  the  number  of  flow  channels 
for  this  set  of  geometries.  The  cross-section  flow  areas  of  both 
geometries  were  0.055  cm  (height)  x  0.09  cm  (width).  Each  dif¬ 
fusion  layer  had  dimensions  of  0.025  cm  (height)  x  14.14  cm 
(width)  x  14.14  cm  (length).  The  number  of  computational  cells 
used  in  the  model  varied  with  complexity  of  the  model.  For 
the  3-channel  flow-field,  the  total  cell  number  was  1.88  mil¬ 
lion  cells,  the  6-channel  had  2.5  million  computational  cells, 
the  13-channel  had  4. 14  million  cells,  the  26-channel  serpentine 
had  6.7  million  cells,  and  the  26-channel  complex  serpentine 
had  5.5  million  cells.  The  large  number  of  computational  cells 
used  in  these  models  could  not  be  solved  using  single  processor 
computer.  Therefore,  PRO-HPC,  a  parallel  computing  imple¬ 
mentation  of  STAR-CD,  with  six  nodes  of  a  Linux  cluster  was 
used  for  these  simulations.  In  the  PRO-HPC  implementation  a 
domain  decomposition  method  was  used  and  each  model  was 
decomposed  into  six  pieces  that  were  sent  to  six  different  proces¬ 
sors  on  the  cluster  for  computation.  Each  processor’s  solutions 
were  communicated  to  the  other  processors  using  MPICH  over 
fast  Ethernet  connections. 

In  this  study,  the  inlet  flow  velocity  was  controlled  by  sto¬ 
ichiometric  numbers  of  1.2  at  anode  and  2.0  at  cathode  and  a 
co-flow  inlet/outlet  configuration  was  used  in  this  study.  The 
operating  pressure  was  101  kPa  absolute  at  the  exit  of  the  cell. 
In  the  model,  the  system  pressure  was  calculated  every  iteration 
to  get  the  correct  inlet  pressure,  velocity,  and  species’  compo¬ 
nent  corresponding  to  inlet  stoichiometry  and  humidification. 
For  temperature  predictions,  the  bi-polar  plate  temperature  was 
maintained  at  70  °C.  The  inlet  humidity  in  this  study  for  both 
the  anode  and  cathode  corresponded  to  80  and  70  °C  dew  point 
temperatures  with  100%  relative  humidity,  respectively.  Note 
that  the  conditions  used  in  this  study  were  the  optimum  PEMFC 
operating  conditions  found  from  internal  experimental  studies. 
All  the  parameters  used  (e.g.,  oxygen  exchange  current  den¬ 
sity,  hydrogen  exchange  current  density,  and  gas  diffusion  media 


Table  1 


Geometry  details,  GDL  properties,  and  MEA  properties 


3 -channel 

6-channel 

13 -channel 

26-channel 

26-channel  complex 

Channel  width  (mm) 

0.9 

0.9 

0.9 

0.9 

0.9 

Channel  height  (mm) 

0.55 

0.55 

0.55 

0.55 

0.55 

Channel  length  (mm) 

3700 

1800 

840 

420 

420 

Rib-spacing  width  (mm) 

0.9 

0.9 

0.9 

0.9 

0.9 

GDL  thickness  (pan) 

250 

250 

250 

250 

250 

GDL  permeability  (m2) 

1.0E-12 

1.0E-12 

1.0E-12 

1.0E-12 

1.0E-12 

GDL  porosity 

0.7 

0.7 

0.7 

0.7 

0.7 

GDL  thermal  conductivity  (W/m  K) 

0.213 

0.213 

0.213 

0.213 

0.213 

GDL  diffusion  adjustment  (%) 

50 

50 

50 

50 

50 

MEA  thickness  (pan)  (including  12.5  p,m  thickness  of  catalyst  layer) 

50 

50 

50 

50 

50 

MEA  thermal  conductivity  (W/m  K) 

0.147 

0.147 

0.147 

0.147 

0.147 

Dry  membrane  density  (g/cm3) 

2.0 

2.0 

2.0 

2.0 

2.0 

Equivalent  weight  of  dry  membrane  (g/mol) 

1100 

1100 

1100 

1100 

1100 

Table  2 

Operating  conditions 
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3-  channel 

6-channel 

13 -channel 

26-channel 

26-channel  complex 

Anode  channel  inlet  conditions 

Gas 

h2 

h2 

h2 

h2 

h2 

Stoichiometry 

1.2 

1.2 

1.2 

1.2 

1.2 

Inlet  temperature  (°C) 

80 

80 

80 

80 

80 

Inlet  relative  humidity  (%) 

100 

100 

100 

100 

100 

Cathode  channel  inlet  conditions 

Gas 

Air 

Air 

Air 

Air 

Air 

Stoichiometry 

2.0 

2.0 

2.0 

2.0 

2.0 

Inlet  temperature  (°C) 

70 

70 

70 

70 

70 

Inlet  relative  humidity  (%) 

100 

100 

100 

100 

100 

Operating  conditions 

Exit  pressure  (kPa) 

101 

101 

101 

101 

101 

Hydrogen  exchange  current  density  (A/cm2) 

0.5 

0.5 

0.5 

0.5 

0.5 

Oxygen  exchange  current  density  (A/cm2) 

0.05 

0.05 

0.05 

0.05 

0.05 

Open  circuit  voltage  (V) 

0.98 

0.98 

0.98 

0.98 

0.98 

Cell  temperature  (°C) 

70 

70 

70 

70 

70 

properties)  were  shown  in  Tables  1  and  2.  In  this  work,  we  con¬ 
sidered  the  effect  of  flow-field  configuration  or  the  path  length  on 
the  PEM  fuel  cell  performance.  The  results  of  local  current  den¬ 
sity,  electrochemical  variables,  temperature,  and  pressure  drop 
with  different  flow-field  configurations  were  compared  and  dis¬ 
cussed  further. 

3.  Results  and  discussion 

Fig.  3  shows  the  polarization  curves  of  the  five  flow-fields  at 
the  conditions  specified  in  Table  2.  The  13-channel  flow-field 
has  the  highest  performance  in  the  IR  region  (i.e.  cell  voltages 
between  0.6  and  0.75)  and  the  mass  transfer  limited  region  (i.e. 
cell  voltages  below  0.6)  followed  in  order  by  the  26-channel 
complex  flow-field,  the  26-channel  flow-field,  the  6-channel 
flow-field,  and  3-channel  flow-field.  In  the  kinetic  region  (i.e. 
cell  voltages  above  0.75),  the  3-channel  flow-field  has  the  best 
performance.  This  high  performance  in  the  kinetic  region  for 


the  3-channel  flow-field  is  related  to  pressure  effects  in  the  flow- 
field  as  will  be  discussed  below.  However,  at  current  densities 
above  1.0  A/cm2,  the  performances  of  both  the  3-channel  and 
6-channel  flow-fields  drop  significantly  compared  to  the  other 
flow-fields.  This  could  be  due  to  the  effect  of  path  length  on 
the  local  performance.  In  order  to  investigate  which  flow-field 
or  path  length  gives  the  optimum  solution,  the  distributions  of 
other  variables  will  be  studied. 

Some  developers  try  to  engineer  their  flow-fields  to  give  uni¬ 
form  membrane  hydration  at  conditions  as  close  to  saturated  as 
possible  without  flooding.  Maintaining  uniform  distribution  of 
membrane  water  content  is  one  method  of  extending  the  lifetime 
of  PEMFC  systems,  because  this  reduces  the  formation  of  local 
hot  spots  and  flooding  that  can  stress  and  damage  ME  As.  Fig.  4 
shows  the  average  membrane  water  content  (A)  at  each  current 
density  for  the  flow-fields  in  this  study.  The  membrane  water 
content  indicates  how  well  membrane  is  hydrated  and  is  the  key 
to  low  membrane  resistance.  Fig.  4  reveals  that  the  3-channel 
flow-field  has  the  highest  membrane  water  content  for  all  cur- 


Fig.  3.  The  performance  curves  of  these  five  flow-filed  configurations  with  the 
condition  of  80  °C/70  °C  dew  point  at  1. 2/2.0  stoich  with  101  kPa  system  pres¬ 
sure  and  70  °C  cell  temperature:  3-channel  (*);  6-channel  (■);  13-channel  (♦); 
26-channel  (A);  26-channel-complex  (x). 
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Fig.  4.  Average  value  of  water  content  of  membrane  at  each  point  on  perfor¬ 
mance  curves  of  these  five  flow-field  configurations:  3-channel  (*);  6-channel 
(■);  13 -channel  (♦);  26-channel  (A);  26-channel-complex  (x). 
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rent  densities,  but  as  current  density  increases  the  membrane 
water  content  decreases.  The  decrease  in  membrane  water  con¬ 
tent  at  increasing  current  density  occurs  due  to  the  increase  in  net 
water  transport  from  anode  to  cathode  with  increasing  current 
density.  The  increased  net  water  transport  from  anode  to  cathode 
results  in  drying  of  the  membrane  because  the  hydration  state  of 
the  membrane  is  a  stronger  function  of  anode  water  activity  than 
cathode  water  activity  [10-12].  Membrane  water  content  has  a 
maximum  of  14.0  for  a  water  vapor  equilibrated  membrane,  but 
for  a  membrane  equilibrated  with  liquid  water  this  value  can 
reach  22.0  [14].  For  the  3-channel  and  6-channel  flow-fields  the 
membrane  water  content  is  above  14.0  for  low  current  densities, 
which  indicates  that  these  flow-fields  have  significant  fraction  of 
liquid  water.  The  presence  of  liquid  water  could  explain  why  the 
overall  performance  of  this  3-  and  6-channel  flow-fields  show 
the  lowest  performance  even  though  they  have  high  values  of 
membrane  water  content.  For  the  flow-fields  with  shorter  path 
length  (i.e.  13-channel,  26-channel,  and  26-channel  complex), 
their  performances  are  consistent  with  the  membrane  water  con¬ 
tent  profile.  The  overall  performance  of  13-channel  flow-field  is 
higher  than  both  of  the  26-channel  flow-fields  because  it  has 
higher  overall  membrane  water  content  at  each  average  current 
density. 

The  water  content  of  the  membrane  is  linearly  related  to  the 
proton  conductivity  of  the  membrane  such  that  higher  membrane 
water  content  leads  to  higher  proton  conductivity.  The  voltage 
loss  caused  by  the  resistance  of  the  membrane  to  proton  trans¬ 
port  can  be  calculated  at  each  point  on  the  polarization  curve  as 
shown  in  Fig.  5.  This  figure  shows  increasing  potential  loss  with 
increasing  of  averaged  current  density.  This  behavior  is  intu¬ 
itive  based  on  the  results  for  membrane  water  content  shown  in 
Fig.  4.  The  most  striking  result  from  this  plot  is  that  the  mem¬ 
brane  in  these  cells  is  only  responsible  for  a  maximum  of  10  mV 
of  the  total  voltage  drop  in  the  cell.  At  each  current  density, 
the  difference  in  membrane  voltage  drop  between  the  differ¬ 
ent  flow-fields  is  only  around  1  mV  even  when  the  membrane 
water  content  is  30%  different.  This  indicates  that  if  flows  are 


Current  density  (mA/cm2) 


Fig.  5.  Ohmic  potential  loss  causing  by  membrane  conductivity  at  each  point 
on  performance  curves  of  these  five  flow-field  configurations:  3 -channel  (*); 
6-channel  (■);  13-channel  (♦);  26-channel  (A);  26-channel-complex  (x). 


well-humidified  then  the  flow-field  should  be  designed  more  for 
eliminating  flooding  conditions  than  trying  to  have  a  perfectly 
uniform  membrane  water  content  distribution. 

The  effect  of  oxygen  concentration  in  the  cathode  on  PEMFC 
polarization  is  reflected  in  the  kinetic  overpotential.  Fig.  6  shows 
the  area  averaged  kinetic  overpotential  as  function  of  current 
density.  This  figure  shows  average  value  of  kinetic  overpo¬ 
tential  at  each  point  on  performance  curves  shown  in  Fig.  3. 
The  3 -channel  flow-field  has  the  smallest  overpotential  at  very 
low  current  densities  because  it  has  the  highest  partial  pres¬ 
sure  of  reacting  gases.  At  low  current  densities  the  amount  of 
water  of  produced  in  the  cathode  is  also  small.  This  combina¬ 
tion  of  effects  leads  to  a  high  partial  pressure  of  oxygen  in  the 
cathode  for  flow-fields  with  longer  path  lengths  even  though 
the  oxygen  stoichiometry  remains  the  same.  This  explains  why 
the  3-channel  flow-field  has  is  the  highest  performance  in  this 
region  as  shown  in  Fig.  3.  For  all  flow-fields,  when  the  aver¬ 
age  current  density  increases  the  kinetic  overpotential  also 
increases.  In  the  IR  and  mass  transfer  regions,  the  13-  and  26- 
channel  flow-fields  have  smaller  kinetic  overpotentials  than  the 
3-  and  6-channel  flow-fields.  When  the  average  current  density 
is  beyond  1000  mA/cm2,  the  value  of  overpotential  is  signifi¬ 
cantly  increased  for  the  3-  and  6-channel  flow-fields  due  to  mass 
transfer.  With  longer  path  length  flow-fields  like  3-channel  and 
6-channel,  the  decreasing  of  gases’  partial  pressure  and  the  accu¬ 
mulation  of  water  are  very  intense  especially  in  the  last  half  of 
the  reaction  area. 

The  pressure  drop  across  the  flow-field  is  an  extremely  impor¬ 
tant  metric  for  PEMFC  system  design  because  it  will  define  the 
types  of  compressors  or  blowers  that  are  needed  and  the  pres¬ 
sure  drop  will  also  have  an  impact  on  the  electrochemistry  in 
the  cell.  Fig.  7  shows  the  effect  of  flow  channel  path  lengths 
on  the  membrane  water  content  and  cathode  pressure  drops  at 
average  current  density  of  1000  mA/cm2.  This  figure  shows  that 
pressure  drop  and  water  content  increase  with  the  increase  of 
path  length.  The  pressure  drop  decreases  as  the  number  of  flow 
channels  increases  because  there  are  more  parallel  paths  that 
the  flow  can  take.  As  the  total  pressure  increases  and  the  mole 


Current  density  (mA/cm2) 

Fig.  6.  Average  value  of  kinetic  overpotential  at  each  point  on  performance 
curves  of  these  five  flow-field  configurations:  3-channel  (*);  6-channel  (■);  13- 
channel  (♦);  26-channel  (A);  26-channel-complex  (x). 
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Fig.  7.  The  effect  of  flow  channel  path  length  from  different  flow-fields  on 
membrane  water  content  and  cathode  pressure  drop  at  average  current  density 
of  lOOOmA/cm2. 

fraction  of  water  remains  constant,  then  the  partial  pressure  of 
water,  and  thus  the  activity,  increases.  The  increasing  activity  of 
water  will  thus  raise  the  membrane  water  content. 

The  distribution  of  reaction  rates  across  the  electrode  surface 
causes  distributions  in  temperature  and  other  electrochemical 
variables  in  a  PEMFC.  Distribution  uniformity  is  important  for 
minimizing  the  material  stresses  on  the  MEA  so  that  its  lifetime 
will  be  improved.  Fig.  8  shows  the  local  temperature  over  the 
anode  membrane  surface  at  lOOOmA/cm2.  The  temperature  of 
the  bipolar  plate  in  all  simulations  is  set  at  70  °C  (343  K)  and 
the  source  term  for  the  heat  of  reaction  is  placed  at  the  cath¬ 
ode/membrane  interface.  The  simulations  show  the  electrode 
temperature  near  the  inlet  region  for  all  cases  is  near  80  °C.  The 


temperature  then  drops  to  around  75  °C  in  the  middle  of  the  flow- 
field  and  maintains  that  temperature  through  the  exit  region.  The 
26-channel  complex  flow-field  has  a  temperature  distribution 
that  is  more  uniform  than  the  other  co-flow  flow-fields  includ¬ 
ing  the  other  26-channel  flow-field.  This  could  be  due  to  the 
special  matching  of  the  maximum  hydrogen  and  oxygen  con¬ 
centrations  at  the  inlet  for  co-flow  configurations.  Furthermore, 
the  average  temperature  between  13 -channel  and  26-channel  is 
only  slightly  different.  This  means  that  there  is  a  limit  to  how 
much  the  temperature  uniformity  can  be  adjusted  by  changing 
the  flow  path  length.  It  is  noted  that  the  temperature  under  the  rib 
areas  is  lower  than  the  adjacent  channel  areas.  This  is  because 
the  heat  produced  in  the  area  under  ribs  is  more  easily  trans¬ 
ferred  to  the  graphite  block  compared  to  heat  produced  in  the 
area  under  the  channels. 

Fig.  9  shows  the  average  temperature  on  the  anode  membrane 
surface  as  a  function  of  the  average  current  density.  The  tempera¬ 
ture  increases  with  increasing  current  density  because  of  higher 
reaction  rate.  As  seen  in  Fig.  8,  there  are  small  differences  in 
average  temperature  at  the  cathode/membrane  interface  for  dif¬ 
ferent  flow-fields,  but  with  the  70  °C  uniform  plate  temperature 
boundary  condition  these  differences  are  only  around  1  °C.  The 
insert  bar  graph  shows  the  temperature  difference  between  the 
maximum  and  minimum  local  temperatures  for  different  flow- 
fields  at  an  average  current  density  of  1000  mA/cm2.  This  graph 
shows  that  longer  path  length  flow-fields  have  larger  tempera¬ 
ture  ranges  than  shorter  path  length  flow-fields.  By  looking  at 
the  temperature  distributions  in  Fig.  8,  it  can  be  observed  that 
the  difference  in  average  temperatures  is  mostly  caused  by  the 
hot  spot  in  the  entrance  region.  The  size  of  the  hot  spot  decreases 
as  the  number  of  channels  increases  and  the  entrance  region  is 
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Fig.  8.  Temperature  (K)  distributions  on  anode  MEA  surface  at  /avg  =  1000mA/cm2:  (A)  3-channel  serpentine  flow-field;  (B)  6-channel  flow-field;  (C)  13-channel 
flow  field;  (D)  26-channel  flow  field;  (E)  26-channel  complex  flow-field. 
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Fig.  9.  Average  value  of  temperature  on  membrane  surface  at  each  point  on 
performance  curves  of  these  five  flow-field  configurations.  Bar  graph  shows 
maximum  temperature  difference  on  membrane  surface  at  7av g  =  1000  mA/cm2: 
3-channel  (*);  6-channel  (■);  13-channel  (♦);  26-channel  (A);  26-channel- 
complex  (x). 

spread  over  more  area.  The  average  temperature  between  13-  and 
26-channel  flow-fields  are  only  slightly  different  which  suggests 
that  there  is  a  limit  to  how  uniform  the  temperature  profile  can 
be  made  by  increasing  the  number  of  flow  channels. 

Fig.  10  shows  current  density  distribution  on  ME  A  surface 
of  the  3-,  6-,  13-,  and  26-channel  serpentine,  and  the  26-channel 
complex  serpentine  flow  paths,  respectively.  The  13-channel 
flow-field  has  a  higher  cell  voltage  than  other  flow-fields.  The  3-, 
6-,  13-,  and  26-channel  flow-fields  give  cell  voltages  of  0.575, 
0.590,  0.620  and  0.605  V,  respectively.  One  important  aspect  of 
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Fig.  1 1 .  Current  density  plots  along  centerline  of  ME  A  from  inlet  toward  outlet 
(i.e.  at  x-direction  at  y  =  7.07  cm)  at  7av g  =  1000  mA/cm2  for  five  different  flow- 
fields. 

these  distributions  is  the  range  of  current  variation.  The  distri¬ 
butions  in  Fig.  10  show  that  the  26-channel  flow-fields  have  the 
least  variation  in  current  density.  The  3-,  6-,  13-,  and  26-channel 
flow-fields  have  differences  between  the  maximum  and  mini¬ 
mum  local  current  densities  of  840,  700,  430,  and  370  mA/cm2, 
respectively.  The  variability  in  the  current  density  can  cause 
local  variations  in  the  cell  temperature,  membrane  conductiv¬ 
ity,  and  water  concentration.  The  non-uniformities  then  cause 
mechanical  stresses  on  the  system  that  reduces  MEA  lifetime. 

Fig.  1 1  shows  further  detail  of  above  contours  on  the  uni¬ 
formity  aspect.  This  figure  presents  the  current  density  plots 
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Fig.  10.  Current  density  (A/cm2)  distributions  at7avg  =  1000 mA/cm2:  (A)  3-channel  serpentine  flow-field;  (B)  6-channel  flow-field;  (C)  13-channel  flow  field;  (D) 
26-channel  flow  field;  (E)  26-channel  complex  flow-field. 
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along  the  centerline  of  ME  A  from  inlet  toward  outlet  (i.e.  at 
y  =  7.07  cm  and  plot  along  x-direction  from  0  to  14. 14  cm.).  The 
plot  of  each  flow-field  shows  regular  peaks  and  valleys  caused 
by  different  rates  of  heat  transfer  in  the  reaction  areas  under  the 
flow  paths  and  under  ribs.  Fig.  1 1  shows  that  the  current  density 
under  the  rib  areas  is  higher  than  under  the  adjacent  channel 
areas.  This  higher  current  density  under  the  ribs  is  caused  by 
local  temperature  as  explained  in  Fig.  8.  Therefore,  the  lower 
temperature  under  the  ribs  compared  to  under  the  channels 
increases  the  membrane  water  content  due  to  the  reduction  of 
saturation  pressure  of  water;  thus,  raising  local  performance.  In 
order  to  reduce  variations  in  current  density  between  channel 
and  rib  areas  the  material  and  thermal  properties  of  materi¬ 
als  including  bipolar  plate,  GDF,  and  MEA  could  be  further 
studied. 

From  above  results,  it  is  noted  that  the  results  between  26- 
channel  conventional  serpentine  and  26-channel  complex  ser¬ 
pentine  are  about  similar  for  all  parameters  and  their  uniformity. 
This  can  be  concluded  that  with  the  same  path  length,  the  perfor¬ 
mance  and  other  variables  will  be  the  same  and  this  has  slightly 
effect  by  flow-field  pattern.  In  these  flow-field  studies,  it  shows 
that  26-channel  give  the  most  uniformity  than  other  flow-fields. 
Furthermore,  the  26-channel  give  the  least  pressure  drop  com¬ 
pare  to  the  other  flow-fields. 

The  simulations  show  that  the  13-channel  gives  the  best  per¬ 
formance  for  a  single  200  cm2  PEMFC.  However,  for  making 
a  PEMFC  stack,  the  26-channel  flow-field  may  be  the  optimal 
choice  to  use  rather  13 -channel  due  to  more  current  density  uni¬ 
formity  and  a  lower  pressure  drop.  Table  3  is  the  summary  of 
these  studies  discussed  above  for  three  points  of  the  V-I  curve 
(i.e.  200,  600,  and  lOOOmA/cm2)  and  it  shows  the  comparison 
of  these  flow-filed  bases  on  channel  length,  cell  voltage,  current 
variation,  temperature  variation,  membrane  water  content,  and 
pressure  drop. 

4.  Conclusion 

The  effect  of  flow-field  path  length  and  pattern  on  the  dis¬ 
tribution  of  current,  temperature,  and  electrochemical  variables 
on  commercial- size  PEMFC  performance  was  studied  using  the 
ES-PEMFC  model.  The  models  used  a  commercial  computa¬ 
tional  fluid  dynamics  solver,  STAR-CD,  and  its  PEMFC  model 
add-on  coded  by  University  of  South  Carolina  for  computation 
in  steady-state,  multi- species,  multi-phase,  and  non-isothermal 
schemes.  The  patterns  studied  were  200  cm2  reaction  area  cells 
with  3-,  6-,  13-  and  26-channel  flow-fields. 

It  is  concluded  that  changing  flow-field  configuration  by  vary¬ 
ing  path  length  of  PEMFC  by  changing  the  number  of  parallel 
channels  can  affect  its  performance  and  uniformity.  The  shorter 
path  length  gives  more  uniform  current  density  distribution  and 
less  condensed  liquid  water  than  the  longer  path.  However,  from 
an  overall  performance  aspect,  the  13 -channel  flow-field  gave 
a  slight  advantage  over  the  26-channel  flow-field  due  mostly 
to  small  differences  in  membrane  hydration.  With  similar  path 
length,  the  performance  appears  to  be  relatively  independent 
of  its  configuration  for  the  26-channel  cases  that  were  tested. 
Therefore,  the  path  length  of  a  PEMFC  flow-field  is  one  of  the 
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major  variables  for  optimizing  the  performance,  efficiency,  and 
durability  of  PEMFC. 
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